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Parametric fluctuations or stochastic signals are introduced into the rectangular pulse sequence to 
investigate the feasibility of random dynamical decoupling. In a large parameter region, we find that the 
out-of-order control pulses work as well as the regular pulses for dynamical decoupling and dissipation 
suppression. Calculations and analysis are enabled by and based on a nonperturbative dynamical decoupling 
approach allowed by an exact quantum-state-diffusion equation. When the average frequency and duration 
of the pulse sequence take proper values, the random control sequence is robust, fault-tolerant, and 
insensitive to pulse strength deviations and interpulse temporal separation in the quasi-periodic sequence. 
This relaxes the operational requirements placed on quantum control devices to a great deal. 

Quantum technology allows one to organize and control the components of a complex system governed by 
the laws of quantum physics'. Control lies at the heart of virtually every aspect of quantum science, and 
in general it is very difficult to achieve. A retrospective and representative study of control follows the 
spin echo effect; inhomogeneous broadening of a spin can be removed by applying a single n inversion pulse 
halfway through the evolution". Quantum control theory has matured over the last two decades, from single spin 
nuclear magnetic resonance experiments to large-scale manipulation during a quantum computation\ As long as 
quantum technology has room to grow, quantum control will remain an active area of research. It has already 
covered a lot of subjects, such as the Lie-algebraic conditions required for strong analytic controllability'' '', 
quantum feedback control on state preparation, stabilization or entanglement improvement' '", time-dependent 
Hamiltonian designation for adiabatic passage"", and exact state transmission'''". Every protocol mentioned 
before needs elaborate design. 

After the investigations based on the control of closed quantum system"* ''^ were gradually extended into the 
open system regime''", active control sequences have been established to defeat the environmental dephasing and 
dissipation effects'""^\ In light of the quantum Zeno effect or similar effects, decoherence of an open quantum 
system could be universally slowed down with ultra-fast modulation"" and concatenated pulses^'. Many of these 
existing works are extremely relying on dynamical decoupling (DD) by subjecting the target system to a series of 
open-loop, high-frequency, control transformations'" ''I DD and its varieties turn out to be a universal and 
efficient class of control method in both theory and experiments. 

The control efficiency of DD, more precisely, the robustness of state fidelity under DD control against 
environment, is mainly determined by the ratio of the characteristic timescale of the environmental correlation 
function and the period or quasi-period of the pulses. In the Markov environment, we are faced with a vanishing 
correlation time. System control is almost impossible since any information contained within a quantum state will 
undertake irreversible loss due to the environment induced error before any external control comes into effect. 
Whereas in a non-Markovian environment, it is possible to drive a quantum state against the environmental 
influence with a properly configured control pulse sequence. A common approximation for DD and its varieties, 
such as UDD (Uhrig DD) and CDD (concatenated DD)'" '"', invokes an ideal delta-function describing a zero- 
width pulse train with unbounded strength. It is essentially a perturbative approach to the effective Hamiltonian. 
In practice, realistic experiments can not be ideal: pulse strength is bounded, duration is finite, and parametric 
fluctuations are inevitable. One can never completely eliminate stochastic quantum fluctuation and envir- 
onmental noise from the laboratory. In other words, the DD mechanisms used in quantum control are generally 
not under total control themselves. In light of this consideration, we are prompted to answer the following 
question: to what extent can a pulse sequence with randomness be effective? 

Traditional DD approaches break down when one discusses the errors randomly occurring in the pulse 
sequence. It is actually the zero-order approximation to a nonperburbative DD pulse sequence'"'. We will instead 
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use a more realistic pulse description in the effective Hamiltonian 
consisting of all of the parameters, such as width, strength, pseudo- 
period and their fluctuations, which is based on the exact quantum- 
state-diffusion (QSD) equation [See Method]. 

Results 

Random Control Model. An illustration of a regularly spaced 
rectangular wave train is shown in Fig. 1(a). This configuration 
allows us to parameterize experimental inaccuracies associated 
with the period t, temporal duration A, and strength O/A. A 
typical example of an imperfect or out-of-order rectangular pulse 
sequence is shown in Fig. 1(b). Notice the time-dependent 
strength, temporal duration, and quasi-period in this case. We will 
use these parameters to measure the fault-tolerance of a system to 
pulse sequence inaccuracies. The regularly spaced sequence 
appearing in Fig. 1(a) should protect the quantum system from 
decoherence and relaxation. The robustness of our scheme can be 
portrayed by its fault-tolerance to fluctuations of three parameters. 
The quasi-period is defined as the time interval between the starting 
point of any pulse and that of its following pulse. The quasi-period 
will generally possess a random distribution. We can express these 
three fluctuating parameters as 



X'=X + Dj(Rand(-l, 1) 



(1) 



where X = t. A, Q>, respectively, D^'s are their individual deviation 
scales and Rand( — 1, 1) denotes a random number uniformly 
distributed between —1 and 1. This is a realistic noise model - 
stochastic fluctuations naturally occur within the quantum 
measurement devices and during the pulse creation process. In the 
limit of a large sample number, M[X'] = X, where M['] denotes the 
statistical ensemble mean. Each random series is instantaneously 
generated and mutually independent. Furthermore, the random 
values are not shared between t. A, and <1> simultaneously. More 
importantly, they are unknown during our random control 
scheme. Without loss of generality, we require A + < t — 
for each rectangular pulse to avoid unnecessary confusion. 

Fault-tolerant Control of Fidelity. Through the derivation by QSD 
equation [See Method]"", a qubit open system with Zeeman split- 
ting CO embedded in a dissipative environment follows an exact 
stochastic Schrodinger equation driven by an effective Hamiltonian: 
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Figure 1 | The diagrams of (a) perfect pulse sequence with constant 
distance, duration, and strength, and (b) irregular pulse sequence with 
random paramters. 



Here c{t) = <1>7A' during the "on" time of random pulse sequence 
which lasts A' and subsequently c{t) = 0 during the "off time of 
duration t' — A' during an interval of instantaneous quasi-period 
t'. Q{t) satisfies an ordinary differential equation dtQ{t) = 
^ + [ — y + ico + ic{t)]Q{t) + Q^{t), and the boundary condition is 
Q(0) = 0. r and y are parameters appearing in the correction func- 

2 

onment coupling strength and l/y is proportional to the 
environmental memory time. When y the non-Markovian 

process is reduced to the Markov limit asymptotically. 

For arbitrary initialization of the system |i//o) = + I'lO); |/'|^ + 
|v|^ = 1, the fidelity that measures the survival probability can be 
expressed by: 



tion of environment G{t,s) 



r is the system-envir- 
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where TZ[-] takes the real part of the following function. One can 
obtain the average of the fidelity over all possible pure state config- 
urations and find 



J^„{t) = - + - 



2jJ&TC|Q(s)] H 

h — 



(4) 



Below we will numerically examine the effects of random control on 
fidelity during time evolution. 

The efficiency of regular and random control is defined in terms of 
the time interval associated with a fidelity drop from unity to 0.95. 
The fidelity is functionally dependent on the bath memory coef- 
ficient y. We illustrate the relationship between the memory function 
and the randomized DD parameters <1), t, and A separately in Figs. 2, 
3, and 4. Each calculation determines the randomization effects 
inherent to a specific control parameter while the remaining para- 
meters remained fixed in a dummy state. The efficiency was then 
evaluated with respect to the ensemble average of these dummy 
states. 

The horizontal axes appearing in Figs. 2, 3, and 4 represent the 
ratio of the stochastic fluctuation scale Dx to the mean value of 
the corresponding parameter X;X=Q>, t, and A, respectively. The 
fluctuations vanish at the origin of each of these coordinate systems. 
This is the condition for regular DD, hence, for a fixed y, the relative 
efficiency of random and regular DD appears in reference to the 




Figure 2 | The moment T when the average fidelity of a single qubit 
reduces to 0.95 from unity as a fiinction of the fluctuation scale for pulse 
strength under different y. We choose <t> = O.lco, t = O.Olojt, A = 0.4t and 
r = CO. 
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Figure 3 | The moment Twhen the average fidelity of a single qubit 
reduces to 0.95 fi-om unity as a fiinction of the period fluctuation scale 
under different /. We choose <l> = 0.2fo, t = 0.02cot, A = 0.4t and T = co. 

intersection of the curve and the vertical axis. These results should 
also be compared to those which naturally emerge in the absence of 
control: T{J^^ =0.95) = lAlcot, 0.&7mt, and 0.65mt when )' = 0.2, 
0.5 and 0.9, respectively. Thus it is evident that both random and 
regular control greatly enhance the survival time of the target state as 
it evolves under open system dynamics; order-of-magnitude 
improvements are found even in the worst situations considered. 

The efficiency has a varying sensitivity to deviations in the DD 
parameters. This behavior is illustrated in Figs. 2, 3, and 4. For 
example, deviations in the pulse strength hardly effect the value of 
T{J-,n = 0.95) . Deviations in the period of the pulse sequence have a 
small impact on the efficiency as well, though a modest decline in 
r(JF„=0.95) tends to occur with increasing D^lx. However, the 
performance of the pulse sequence of DD becomes sensitive to devia- 
tions in A as the bath approaches the strong non-Markovian regime. 
This sensitivity is exemplified in Fig. 4 by the 40% reduction of Tover 
the evaluation scale 5A/4 < t' < 15A/4 with y = 0.2, a typical strong 
non-Markovian condition. These results indicate the effectiveness of 
noisy quantum DD. The relative performance of regular and random 
control depends on the environmental coefficient y, improving as the 
environment becomes more Markovian. For near-Markovian pro- 
cesses, little distinction can be made between completely regular and 
completely random control. 

The common thread among these three figures is the memory 
effect of the environment. In terms of the fidelity expression in Eq. 
(4), strong suppression of dissipation requires the exponential func- 



tion Q[t) =exp 
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to approach unity in the presence of 



c{t). In absence of control pulse, 
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where y = y — iw. A straightforward derivation shows that when y is 
dominant in Eq (5), the absolute value of Q(t) approximately decays 
in an exponential way with time. Thus the control pulse is not able to 
dynamically decouple the open system for dominantly large y, i.e. 
within the near-Markovian environment. On the contrary, smaller y 
implies a longer memory and a slower characteristic decay of the 
quantum state. Coherence can be preserved in this case with a control 
pulse of proper frequency. Appropriate pulse sequences could 
roughly wash out the dissipation influence. In summary, we have 
found that standard quantum DD can be quite effective when fluc- 
tuations occur in the pulse sequence. The performance is largely 
insensitive to fluctuations in pulse strength and sequence period. It 
is, however, quite sensitive to several other factors. The greatest 



Figure 4 | The moment T when the average fldeUty of a single qubit 
reduces to 0.95 from unity as a function of the temporal duration 
fluctuation scale under different y. We choose <1) = 0.2co, t = 0.02ojf, A = 
0.4t and F = co. 

influence, to no surprise, comes from the structure of the envir- 
onmental noise process itself A non-Markovian environment is 
required for effective DD in both cases, regular and random. 
Secondly, the control pulse sequence requires a sufficiently fast pulse 
series, i.e. a sufficiently small average value of t. Dissipation effects 
from the system-environment interaction create disordered system 
dynamics, these effects can be neutralized if a sufficiently large num- 
ber of control photons interact with the system during the character- 
istic correlation time of the environment. Hence, the required pulse 
rate depends on the fixed bath correlation time, a time that is inver- 
sely proportional to y. 

A third condition, which has not been discussed thus far, concerns 
the ratio between the pulse width and pulse periodicity. This depend- 
ence is clearly illustrated in Fig. 5 where comparisons in the fidelity 
evolution are made between regular and stochastic control for vari- 
ous A/t. Perhaps the most striking comparison results when A/t = 
0.3 - random control can actually outperform regular control in 
some occasions. To be fair, we concede that regular DD is more 
favorable towards the beginning of this example evolution. As the 
interpulse separation diminishes, e.g. A = 0.75t, it becomes difficult 
to discriminate the curves and the fidelity remains roughly the same 
in either case. Moreover, in Fig. 6, similar patterns emerge when we 
examine inter-parameter dependencies between A and t. Both of 
these figures correspond to the same average A and t. Again, we find 
examples where random control permits longer survival times than 




Figure 5 | The time evolution of the average fldehty of a single qubit 
under regular control (blue hnes) and under random control (red hnes) 
with different A. We choose y = 0.3, <S> = 0.2co, t = 0.02cot and T = co. 
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the regular case for long evolutionary periods. However, in this case 
random control outperforms regular DD control in early evolution- 
ary intervals as well. Over the entire evaluation interval, random 
control yields a fidelity greater than 0.9 when the ratio between 
the average values of A and z is optimized (A/t S 0.4). Although 
the optimized random controls yield slightly lower fidelities than the 
regular controls for these average values, they could adequately 
approximate the ideal regular sequence if the performance require- 
ments are not prohibitively restrictive. 

Discussion 

The effectiveness of the random DD pulse sequence is fairly state- 
independent. The similarity in performance for several initial qubit 
configurations is illustrated in Fig. 7. These configurations are 
defined in terms of the initial population of the upper level of 
the open qubit system. We plot the fidelity as a function of |^|^ rather 
than ji in light of Eq. (3). Curves associated with a particular |/(|^ 
represent the common result of a group of quantum states with the 
same upper level population. We find a slight decrease in fidelity 
control efficiency with increasing Nevertheless, random control 
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Figure 7 | The time evolution of the fidelity of a single qubit under 
random control witli different initial states indicated by population \//\ ^. 
We choose y = 0.3, <S> = 0.2co, t = 0.02ft)t, A = 0.4t, Da = 0.2t, = 0.2t 
and r = CO. 



provides robust and reliable stabilization against damping effects for 
all possible initializations. In the short time regime, the efficiency of 
random control is nearly the same for every possible initial pure state. 

We have introduced a general unknown stochastic disturbance 
into the nonperturbative quantum DD model to account for imper- 
fect laboratory control. An error-free control sequence consists of a 
sequence of identical, equally spaced pulses. In our nonperturbative 
method, this sequence has already far beyond the traditional DD 
pulse: the bang-bang (BB) control sequence. With high precision, 
we can control the average values of the parameters which character- 
ize the sequence, but we cannot control the errors which occur on 
them or predict any occurrence in advance. Using the quantum- 
state-diffusion equation, we have analyzed the performance of a 
realistic quantum DD sequence prepared with random laboratory 
parametric fluctuations, and find it to be effective in controlling the 
dynamics of an open quantum system over a wide range of non- 
Markovian processes. 

The control efficiency, which is measured by the average fidelity 
over initial pure states, is found to be dependent on the average 
frequency of the control pulse, the environmental memory, and 
the ratio of the average pulse duration time A to the period t. 
When the mean values of A and t are optimized into a very loose 
condition, random dynamical decoupling sequences can adequately 
approximate the ideal regular sequence if the performance require- 
ments are not prohibitively restrictive. This will alleviate the experi- 
mental requirements placed on pulse sequence generation in 
quantum control to a great extent. 

The survival time, which is measured by the time moment T when the 
system fidelity reduces to a given value J- from unity, is associated with 
the quantum speed limit (QSL) of open quantum systems* The lower 

bound of the required time of evolution is T = 2(cos^'jF)^y' ('^^^), 

where S = t-' 0)|.//o>('/'ol^I(f. 0)||"- T^-^' are the 
Kraus operators characterizing the time evolution of the open 
system, and ||j4|| is the Hilbert-Schmidt norm of A. For our dis- 
sipation model, the summation in the expression of S is 

|/.nQW|||Re[Q(0]||Q(OI + |Q(f)lVl/l/'l'-l + lQ(Ol'}. where 

the larger the population of the upper state is, the shorter the evolu- 
tion time T is. It can also be expected that if a system under a certain 
control yields nearly vanishing Q{t) during the evolution, t could be 
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very large such that T is greatly enhanced by the control. Although 
our system is in a non-Markovian environment, which is beyond the 
scope of Lindblad master equation, yet it turns out that when j/ij^ = 

0. 5 and there is no external control applied to the system, the survival 
time r corresponding to .F = 0.95 is around 0.13 (in the same unit as 
in Fig. 7), whereas the lower bound t ~ 0.05 which serves as a 
significant lower bound for T. In general, the QSL under non- 
Markovian dynamics and external dynamical control remains an 
interesting topic for future research. 

The randomness considered here is entirely different from the 
random decoupling schemes initiated by L. Viola et al. (see 
Refs.'""'*''). Based on standard (ideal) n pulse, the random DD pro- 
posed in these works allows the "control propagator" to follow a 
random but known path on the group of rotating operations. Both 
the past "control operations" and the moments at which they are 
applied are known, but the future "control path" is random. It should 
not be confused with the stochastic error considered in the phase or 
rotation axis of n pulse of DD, since by which the connection 
between the control efficiency and the randomly errors in shape of 
pulse is still not clear. Our DD with random control assumes the 
occurrence of unknown stochastic fluctuations over pulse sequence 
parameters, including duration time, period, and strength. The mag- 
nitude of the errors or fluctuations are so big that the previous DD 
pulse sequence now has not much significance than a reference. This 
investigation also differs from the control with no control procedure 
proposed by some of the authors of this work, see Ref 46. There, an 
out-of-control environment consisting of disordered external white 
noise was shown to suppress decoherence in a more ordered system. 
Here, we demonstrate the effectiveness of quantum state stabilization 
when irregular control pulses are applied. 

Methods 

We employ an exact stochastic Schrodinger equation, the QSD equation, to invest- 
igate this stochastic special DD process. This allows one to include the function of 
external pulse sequence directly into the microscopic quantum model. For an arbit- 
rary model (setting h ~ 1): 

-HLot=Hsys(t)+ {Sk^'^k + '^-^-) + (6) 

k k 

where a], and ak respectively denote the creation and annihilation operators for the k- 
th mode in the bosonic bath. In the interaction picture with respect to ^^k ^^^^i^^' 
exact QSD equation generally reads'*^ "'**; 

d,il,,{z') = -i[H,y,{t) + iLz; -iL*d{t, z*)]iA,(z*) = -iHetf(Oi/',(z*)- (7) 

Here Hsy^it) is the system Hamiltonian that could absorb an arbitrary DD pulse 
function. We have dropped the ket notation for = | L is the coupling 

operator of the system with the environment. Zj = — i gl zl e""'' ' is the envir- 
onmental Gaussion noise process satisfying M[z*] — M[z*2*] —0 and 
M[z[2*] = G{t, s) and G(f, s) is the environmental correlation function. 0{t, z*) 
includes system operators and the environmental noises 

0(f, z*)i//f(z*) — ii5G(f, 5)0(f, 5, z*)t//,(z'') with O-operator satisfying; 
Jo 

d,0{t,s,z') = [-iHMt),0{t,s,z')]-L^'^^^, (8) 

and 0(5, s, z*) = L. The states of the system correspond to a particular 

"configuration" z*, thus the system density matrix is recovered by Pf — M[|i//f(2*)> 

Consider the random control over a single two-level system (qubit) which is under 
the dissipative influence from the environment, so that H^ys — E{i)oJl and L — o— 
The non-perturbative pulse sequence indicated by a function of time c{t), regular or 
random, contributes to the system energy shift E{t) — co + c{t). For this model, the 
exact effective HamUtonian in Eq. (7) is found to be Eq. (2). 
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